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ABSTRACT 

We have detected circumstellar absorption lines of the ^^CN and ^^CN Violet and 
Red System in the spectrum of the post-AGB star HD 56126. From a synthetic spec- 
trum analysis, we derive a Doppler broadening parameter of 6 = 0.51 ± 0.04 km s~^, 
^^CN/^^QN^ 38^2, and a lower limit of 2000 on ^^CN/^^CN and ^^ci^N/^^qis^^ ^ gj^, 
pie chemical model has been computed of the circumstellar shell surrounding HD 56126 
that takes into account the gas-phase ion-molecule reaction between CN and C"*". From 
this we infer that this reaction leads to isotopic fractionation of CN. Taking into account 
the isotopic exchange reaction and the observed ^^CN/^^CN we find ^^C/^^C^^ 67 (for 
2kin = 25 K). Our analysis suggests that ^^CN has a somewhat higher rotational tem- 
perature than ^'^CN: Trot = 11-5 ± 0.6 and 8.0 ± 0.6 K respectively. We identify possible 
causes for this difference in excitation temperature, among which the A^" dependence of 
the isotopic exchange reaction. 



Subject headings: line: identification - molecular data - molecular processes - 
stars: AGB and post-AGB - stars: circumstellar matter- individual stars: HD 
56126 



1. Introduction 

The evolution of a low-mass red giant terminates 
on the tip of the Asymptotic Giant Branch (AGB). 
The remnant, which is a carbon-oxygen core with a 
dilute extended convective envelope, evolves along a 
constant luminosity track toward the White Dwarf 
(WD) phase. This transition from the AGB to the 
WD is referred to as the post-AGB phase, or the 
Pre-Plancrtary Nebulae (PPN) phase. Since low-mass 
stars are the main contributors of carbon, nitrogen, 
and s-process elements to the interstellar medium 
(ISM) (Forestini & Charbonnel 1997), accurate in- 
formation of the chemical composition of AGB and 
post-AGB stars is important. Post-AGB stars may 
provide useful information on the composition of the 
gas returned to the ISM. In particular, the fact that 
their spectra are much simpler than those of the cool 
AGB stars affords novel opportunities to infer details 
of the composition of AGB stellar envelopes. 

The present photosphere of HD 56126 {T^s = 7000 

K) is carbon-rich (C/0~1.4), metal-poor ([Fe/H]~- 
1.0), and enhanced in s-process elements ([s/Fe]=1.7) 
(Klochkova 1995). Such an abundance pattern re- 
sembles that of carbon stars (Lambert et al. 1986, 
Utsumi 1970) that result from the dredge-up of nu- 
cleosynthesis products as a star evolves up the AGB. 
HD 56126 is the prototype of a group of post-AGB 
and AGB stars that exhibit absorption lines from cir- 
cumstellar C2 and CN (from now on C means ^^C 
and N means -'^^N) in their spectra (Bakker et al. 
1996 (Paper I), Bakker et al. 1997 (Paper II)). In or- 
der to facilitate further discussion, we introduce the 
term "C2CN" stars to refer to this group of stars. 
Molecules are present in the detached shell surround- 
ing the star and the chemical composition of the dusty 
shell reflects the photospheric composition of the star 
when it was at the tip of the Asymptotic Giant Branch 
(TP- AGB). The molecular composition of the shell 
has evolved, for example a simple molecule like CN 
is believed to be formed as a result of photodissocia- 
tion of complex molecules (in this case HCN) by the 
interstellar radiation field. 

Excitation of the CN molecule's X^S+ ground 
state is likely controlled by absorption and re-emission 
of photons in the Violet (B'^'E+ X'^'£+) and Red 
(A^n X^S+) System, and pure rotational tran- 
sitions (principally de-excitation) in the X^S+ state. 
At the low densities of the circumstcUar shell, the pure 
rotational transitions will cool the rotational ladder 



below the gas kinetic temperature (sub- thermal). For 
a symmetric molecule without a dipole moment (e.g. 
C2), pure rotational transitions are forbidden and the 
molecule can not efficiently cool. The excitation tem- 
perature is higher than the kinetic temperature of the 
gas (supra-thermal). Through measurement of the 
level population and modeling of the (de-)excitation 
processes, the CN molecule serves as a probe of the 
physical conditions in the shell. A molecule like CN 
additionally affords an opportunity to compare col- 
umn densities of isotopomers, e.g. the ratio of CN 
to ^^CN. In turn, understanding of the chemistry en- 
ables the isotopic ratio ^^C/^^C to be derived from 
the measured ^^CN/^^CN ratio. 

In this paper (Paper III) we present and discuss 
the first measurement of the CN/^'^CN ratio and an 
estimate of the ^^C/^^C ratio and lower limits on 
i^C/i^'C and "N/^^N in the circumstellar shell of a 
post-AGB star. In Sec. 2 we discuss the observa- 
tions and the dataset of equivalent widths used in our 
analysis. Sec. 3 goes into the details of the molecu- 
lar parameters used, and Sec. 4 describes the analysis 
and results. A discussion of the results is presented 
in Sec. 5. 

2. Observations and equivalent widths 

High- resolution spectra were obtained with 2.7 m 
Harlan J. Smith telescope of the W. J. McDonald ob- 
servatory and the coude cross- dispersed echelle spec- 
trograph (Tull et al. 1995). Spectra at a resolu- 
tion R = A/AA ~ 140,000 were acquired of the 
(w' - v") = (0 - 0) band of the CN Violet System and 
R ~ 200, 000 of the (3-0) and (4-0) bands of the CN 
Red System. Observations were made from December 
1996 to March 1997 (Table 1). The spectral resolu- 
tion of the spectra was determined from the FWHM 
of the emission lines in the accompanying ThAr arc 
spectrum and corrected for the intrinsic width of the 
ThAr emission lines. 

Observations of 30 minutes have been co-added to 
give the final spectra for each night. All spectra were 
reduced using IRAF in the standard manner. Individ- 
ual integrations were corrected for the velocity shift 
due to the Earth rotation and instrumental drift be- 
fore combining them to an average spectrum. The 
final S/N ratio of the spectra is 35 for the Violet (0- 
0) band and 200 and 90 for the Red (3-0) and (4-0) 
bands respectively. 

Portion of the final spectra are shown in Figs. 1 
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Fig. 1. — CN Violet System (0,0) band towards HD 56126. The upper panel shows the observed spectrum and 
a rectified spectrum corrected for the underlying photospheric features. The rectified spectrum contains only 
circumstellar lines. The lower panel shows on an expanded wavelength scale the strongest line (the Ri, R2, and 
^Q2i blend for N" = 1), and demonstrates that we have not detected ^^CN and C^^N. The dashed spectrum is 
a synthetic spectrum computed for b = 0.51 km s^^ , Trot = 11-5 K for CN, T^ot = 8.0 for the CN isotopes, and 
the isotope or lower limit isotope ratios as determined in this work. The solid spectrum is the synthetic spectrum 
convolved to our spectral resolution of i? = 140, 000. 
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Fig. 2. — CN Red System (4,0) band towards HD 56126. The upper panel shows the observed spectrum (offset 
by 0.2), a rectified spectrum corrected for the underlying photospheric features, and convolved synthetic spectrum 
(offset by -0.2) The rectified spectrum contains only circumstellar lines. The lower panel shows on an expanded 
wavelength scale centered around the strongest line (the ^Q2i(l) or R2(l) blend) (lower panel). The dashed and 
solid synthetic spectrum are explained in the caption of Fig. 1. The ^^CN (4,0) band is red-shifted by about 54 
A and not covered by this spectrum. The solid spectrum is the synthetic spectrum convolved to our spectral 
resolution of i? = 200, 000. 
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Table 1: Log of observations. 



Date 


HJD 


Int. [s] 




Ac [A] 


SNR 


Remark 


29 Dec. 96 


2450446.7158 


6x1800 


200,000 


7558 


100 


Red (3,0) + Ha 


02 Jan. 97 


2450450.8708 


3x1800 


140,000 


4062 


20 


Violet (0,0) + Ca II H + He 


28 Feb. 97 


2450507.5840 


6x1800 


140,000 


4056 


15 


Violet (0,0) + m 


01 Mar. 97 


2450508.5696 


6x1800 


140,000 


3874 


22 


Violet (0,0) + H(5 


03 Mar. 97 


2450510.6087 


8x1800 


200,000 


6086 


70 


Red (3,0) + Red (4,0) 


Co-added spectra 












Violet System (0,0) 


7.5 hrs 


140,000 


3877 


35 




Red System (3,0) 


7.0 hrs 


200,000 


6945 


200 




Red System (4,0) 


4.0 hrs 


200,000 


6195 


90 





and 2. Prom each final spectrum we removed the 
molecular features and fitted a high order spline that 
followed all photospheric features. The observed spec- 
trum what divided by this fit to obtain the recti- 
fied spectrum which only contains the molecular fea- 
tures. Circumstellar CN lines are readily distinguish- 
able from the much broader photospheric lines of var- 
ious atomic species. The equivalent widths (Tables 2 
and 3) of the molecular lines were measured relative 
to the local continuum (taking into account stellar 
absorption lines). A search for ^^CN and C-'^^N of the 
Violet System (0,0) band was unsuccessful: an upper 
limit to their equivalent width is 1.5 mA. 

HD 56126 is a pulsating star. The photospheric 
spectrum from the star might therefore change from 
one run to the other and combining the spectrum 
from different runs could lead to errors. We there- 
fore have measured the equivalent width of the lines 
for each night separately and averaged to obtain the 
final equivalent width. The spectra in Figs. 1 and 2 
are combined using all available data and there is no 
indication that the photospheric spectrum has indeed 
changed from one run to the next. 

This high quality dataset was extended with pre- 
vious data of the Red System (1,0), (2,0), and (3,0) 
bands (presented in Paper I), provided by spectra col- 
lected with the WHT/UES at a spectral resolution of 
R K. 50, 000. Our new measurements for the Red (3- 
0) band are in good agreement with measurements of 
Paper I. The weaker lines were only detected at our 
higher resolution. Our new measurements for the Red 
(3-0) band are in good agreement with measurements 
of Paper I. 



3. Molecular data 
3.1. CN Violet System 

Wavelengths for CN Violet System (0,0) were com- 
puted from the wavenumbers of Prasad et al. (1992) 
and the isotopic wavelengths for -'^^CN Violet System 
(0,0) from Jenkins & Wooldridge (1938). The iso- 
topic shift of the ^^CN Pi(iV" = 1) blend was ob- 
tained from our observed spectrum. Lines positions 
not listed in these papers were computed using the 
molecular constants given by Prasad et al. and the 
standard relations for the mass dependence of the 
various molecular constants (cf Bernath 1995). The 
computed wavenumbers were in good agreement with 
the laboratory measurements. Predictions were ex- 
tended to the isotopes ^'^CN and C^^N. Table 2 lists 
the wavelength shift of the various isotopes relative to 
CN. The isotopic shift is the same for the three lines 
within one blend. For example: the ^^CN ^Q2i (2) 
line is at 3873.593 A. In all cases the index of refrac- 
tion of standard air was computed using the formula 
given by Morton (1991). 

The /v-value of the Violet System is well de- 
termined experimentally and theoretically. Theo- 
retical calculations by Knowles et al. (1988) and 
Bauschlicher, Langhoff & Taylor (1988) predict /v(0- 
0)=0.0345 and 0.0335 respectively. An earlier cal- 
culation by Larsson, Siegbahn & Agren (1983) gave 
/v (0-0) =0.0324. The most accurate experimental re- 
sults are the measurements of the radiative lifetime of 
the B^I]+ state. Low rotational states of the v' = Q 
level of the B state have measured lifetimes (in ns) of 
60.8 ± 2.0 (Luk & Bersohn 1973), 65.1 ± 0.8 (Jack- 
son 1974), and 66.7 ± 1.4 (Duric, Erman & Larsson 
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Table 2: CN Violet System (0,0) with /(o,o) = 0.0330. Derived column densities are given in Table 4. 
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0.3481 
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0178iVf6') 
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< 2.0 
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OO / U.\JU 1 i 
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0009 










Ri(5) 
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0.3001 


0.550c 


0.218c 
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< 2.0 
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0.1921 
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0.137c 
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3874.6021 


0.1761 
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0.119c 
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0.0110 










Pifi") 


3875.7591 


0.130o 
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0.083c 


0.0110 
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3876.830c 








0.0009 










Pi(3) 


3876.8341 


0.0631 


0.130c 


0.050c 


0.0141 


0.0141iV(3) 


27.6 ±0.3 


4.2 ±0.5 




P2(3) 


3876.8361 








0.0132 










''Ql2(4) 


3877.336c 








0.0006 










Pi (4) 


3877.3411 


0.0561 


0.091c 


0.035c 


0.0147 


0.0147iV(4) 


24.1 ±0.3 


< 1.0 




P2(4) 


3877.3451 








0.0141 










^Ql2(5) 


3877.821c 








0.0003 










Pi(5) 


3877.8281 


0.028c 


0.054c 


0.020c 


0.0150 


0.0150iV(5) 


3.9 ±0.5 


polluted 


perturbed 


P2(5) 


3877.8251 








0.0147 












3877.8851 












< 1.0 




perturbation 


^Ql2(6) 


3878.286c 








0.0002 










Pi(6) 


3878.2941 


0.008c 


0.018c 


0.006c 


0.0152 


0.01527V(6) 


< 1.0 


< 1.0 




P2(6) 


3878.3001 








0.0150 











a: isotopic shift relative to CN. Example: R2(l) "CN at 3873.996+0.348=3874.344 A 



c: computed; 1: laboratory data; o: spectrum 
bl: blended with the ^^CN P(5) blend 

For lines with laboratory data: v(CN)= 77.5 ± 0.5 km s'^ and t;(i3cN)= 76.7 ± 1.8 km s'^ 
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Table 3: CN Red System (4,0) with /(4.0) = 1.09 x 10"^. Derived column densities arc' giv(m in Table 4. 

B{N") A,,,i, [A; /a-'j-.a"./" X 10^ ./;;ir X l(.)^.V(A''0[cm--] Ua [mA] Rom^~ 



'^RaiCS) 6186.2931 0.101 0.058N(3) < l.S' 



"R2l(2) 


6187.7641 


0.112 


0.067N(2) 


1.8 ±0.5 


"R2l(l) 


6189.3971 


0.127 


0.085N(1) 


3.1 ±0.5 


^R2l(0) 


6191.159c 


0.151 


0.151N(0) 


2.4 ±0.1 


R2(3) 


6192.0871 


0.202 






^Q2l(3) 


6192.101c 


0.217 


0.211N(3) 


< 1.5 


R2{2) 


6192.2551 


0.203 






^Q2l(2) 


6192.274c 


0.237 


0.223N(2) 


5.3 ±0.6 


R2(l) 


6192.572c 


0.215 






^Q2l(l) 


6192.6051 


0.264 


0.248N(1) 


8.0 ±0.2 


^Q2l(0) 


6193.0081 


0.363 


0.363N(0) 


7.2 ±0.1 


Q2(l) 


6194.4851 


0.363 


0.182N(1) 


6.7 ±0.1 


QP2l(l) 


6194.550c 


0.091 






Q2(2) 


6195.4671 


0.316 


0.188N(2) 


4.5 ±0.5 
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6196.6221 


0.323 


0.197N(3) 


1.7 ±0.2 


QP2i(3) 


6196.638c 


0.103 







V2C2) 6197.4.381 0.091 0.03GX(2) < 1.5 



1978). Theoretical values (including the small con- 
tribution from the B^S+ ^A^II) arc 60.7 (Knowles 
et al.), 62.4 (Bauschhcher et al.) and 66.6 (Larsson 
et al.) We adopt /v(0 — 0) = 0.033 which is approx- 
imately the average of these various theoretical and 
experimental estimates. The uncertainty is not more 
than a few per cent. 

3.2. CN Red System 

Accurate wavelengths of Red System CN lines were 
taken from Davis & Phillips (1963), and the iso- 
topic shifts for ^'^CN Red System from Hosinsky et 
al. (1981). Wavelengths not found in literature were 
computed using the molecular constants of Prasad & 
Bcrnath (1992) and compared to those of the SCAN 
tape (J0rgensen & Larsson 1990). 

Theory and experiment have not yet fully con- 
verged for the Red System /p-vahic - sec brief re- 
view by Larsson (1994). Theoretical predictions are 
/r(3 - 0) = 3.34 X 10-^ (Knowles et al. 1988), 
3.35x10"'' (Bauschlicher et al. 1988), and 4.58 x lO"'' 
(Larsson et al. 1983). Davis et al. (1986) measured 
the /r- values of six Red System bands relative to the 
Violet (0-0) band using absorption lines produced by 
the transmission of light through a column of CN 
in a furnace. Their result, adjusted by 3% to ac- 
count for the difference between their adopted and the 



above recommendation for /v(O-O), gives /r(3 — 0) = 
(2.9 ± 0.2) X 10~^, a value slightly smaller than the 
most recent theoretical calculations. Davis et al.'s re- 
sults for other bands were similarly smaller than the 
theoretical predictions. Gredel, van Dishoeck & Black 
(1991) analyze observations of CN Violet (0-0) and 
Red (2-0) interstellar absorption lines in the same line 
of sight using the /p-valucs of /v(0 — 0) = 0.0342 and 
/r(2 - 0) from Davis et al. (1986). Gredel et al. re- 
mark that the Violet and Red lines give the same col- 
umn density but "a slightly higher Red System oscil- 
lator strength would improve the agreement between 
the Violet and the Red System data" . Inspection of 
Gredel et al.'s results suggests that the theoretical 
value given by Knowles et al. and Bauschlicher et al. 
fit the suggestion of "a slightly higher" /^-value but 
Larsson et al.'s value is probably too large. There 
remains an unresolved question posed by measure- 
ments of the radiative lifetimes of vibrational levels 
of A^n state. What appear to be the best measure- 
ments, radiative lifetimes from laser-induced fluores- 
cence (Lu, Huang & Halpern 1992), are consistently 
shorter than the theoretical predictions: the differ- 
ence amounts to about 30 % for v'=0 but a factor 
3 for v'=7. These lifetimes imply that our adopted 
/ri-values may be underestimates. We adopt the /r,- 
value given by Knowles et al. and Bauschlicher et 
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al. with a small correction for our adopted fy value: 
/r(0 - 0) = 23.7 X 10-^ /r(1 - 0) = 19.1 x lO"", 
/r(2 - 0) = 9.0 X 10-^, /r(3 - 0) = 3.3 x 10"*, and 
/r(4-0) = 1.1x10-4. 

3.3. Oscillator strength of individual lines 

The oscillator strength of a transition is given to an 

acceptable precision by /a?/ J/ jv" J" ~ {vn' j' .N" j" /vq) 
fv'v"Sj'j"/ (2 J" + 1) where iyN'j',N"j" and i/q are 
the frequency of the transition and the band origin 
respectively, fv'v" is the band oscillator strength 
discussed in Sect 3.1 & 3.2, and S'j'j" is the Honl- 
London factor (also called the natural line strength). 
The computation of the Honl-London factors is de- 
scribed in App. A. 

Almost all CN "lines" consist of two or three un- 
resolved transitions. For example, the R-branch lines 
of the violet system (Table 2 and Fig. 1) are effec- 
tively a triplet of transitions (Ri(A^"), R2(iV"), and 
^Q2i(A''")). For the Red System, several branches 
provide an accompanying satellite line (Table 3 and 
Fig. 2), e.g. , R2(A^") is accompanied by ^^Qsi (TV"). 

The determination of column densities may use a 
curve of growth computed for a single line. If unre- 
solved components of a blend are separated by more 
than about the Doppler width AAd, the equivalent 
width of the lines of individual transitions are added. 
If the separation is much less than AAd, the blend 
may be presented as a single line with a so-called ef- 
fective oscillator strength (see App. A). For separa- 
tions comparable to AAd, there is no simple way to 
represent the blend on the curve of growth for a single 
line. 

A correct treatment of the overlap of spectral lines 
in a blend necessarily requires the computation of syn- 
thetic spectra. A computer program (MOLLEY-CN) 
was developed that computes wavenumbers, Honl- 
London factors and oscillator strength for any (iso- 
topic) band of the CN Red and Violet System. From 
these numbers a synthetic spectrum based on a Voigt 
profile is computed. Input to the program is the rota- 
tional temperature, total column density, isotope ra- 
tios, Doppler 6-parameter, and the spectral resolution 
of the output spectrum. MOLLEY-CN is available on 
request from the authors. As an aid to those who pre- 
fer the method of effective oscillator strength, we have 
investigated the error this latter method introduces 
on the derived column densities. For this purpose 
we have chosen the CN Violet (0,0) R(2) blend for 



which the two principal transitions are separated by 
0.54 km s"^. For a range of Doppler 6-parameters the 
column density of the N" = 2 (Fi and F2) level has 
been computed using the effective oscillator strength 
and synthetic spectra for a line with an equivalent 
width of 37.1 mA. Fig. 3 demonstrates the results. 
The derived 5-parameter in this work is6=:0.51±0.04 
km s"^, and the effective oscillator strength method 
is clearly unacceptable for this work. We therefore 
will use the method of computing synthetic spectra 
to determine the column density responsible for the 
observed equivalent width of a "blend" . The only re- 
maining limitation is that the line might be blended 
with transitions originating from other N" level and 
that the rotational temperature must be given to com- 
pute such a blend. This limitation is easily overcome 
by iteration using an improved estimate of the rota- 
tional temperature. 

Fig. 3 can be qualitatively understood in the fol- 
lowing way. Use of fcs assumes that the lines which 
form a blend (R(2) is a triplet) have exactly the same 
wavelength. This approximation is valid only for b- 
values much larger than the separation between the 
lines of a blend. Fig. 3 demonstrates this as the com- 
puted column density from the two methods are the 
same for large 6-values. For small 6-values, the use of 
/off fails because the lines are (partly) separated and 
the blend absorbs more efficiently than assumed by 
the fee method. For a given equivalent width of a 
line, this means that the spectrum synthesis method 
yields a lower column density than the /off method. 
The R(2) blend has three lines. For moderate column 
densities, the contribution of the weakest line ('^'Q2i 
(2)) can be neglected and the ratio of the column 
density as derived from the two method is given by 
the inverse of the slope of the curve of growth. The 
slopes are: 1 for the linear part, 0.5 for the damping 
wings, and < 0.5 for the saturated part of the curve 
of growth, and the maximum ratios are therefore: 1, 
2 and > 2. We see that in Fig. 3 the ratio goes to 
2 when the two strongest lines of the triplet are on 
the saturated part of the curve of growth. For larger 
column densities, the blend is highly saturated and 
the broad damping wings contribute significantly to 
the equivalent width. The damping wings will be the 
same in both cases and the ratio will go to one. We 
note that in the "real" world lines never become this 
strong. 
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4. Analysis 

4.1. Velocities 

Prom the molecular lines with tabulated labora- 
tory wavenumbers, we find an average heliocentric 
radial velocity of 77.5 ± 0.5 km s^^. To within the 
errors of the measurements, the velocity is indepen- 
dent of the rotational level A^". All velocities given 
in this article are corrected to the heliocentric rest 
frame. HD 56126 is a pulsating star (Oudmaijcr & 
Bakker 1994, Lebre et al. 1996, discussion of Pa- 
per 11) and the photospheric velocity derived from 
our spectra will not be a good estimate of the sys- 
temic velocity of the star. In order to derive an ex- 
pansion velocity wc use the systemic velocity derived 
from CO radio line emission (for overview see Paper 
I). For w*,co = 85.6 ± 0.5 km s^^ wc arrive at an 
expansion velocity of Woxp = 8.1 ± 0.7 km s~^. This 
result is consistent with that derived from the width 
of the CO lines: Vexp = 12.1 ± 1.0 km s~^ (Nyman 
et al. 1992). The expansion velocity is similar to ex- 
pansion velocities of AGB stars suggesting that the 
molecule-containing layer around the post- AGB star, 
which is possibly part of the terminal "superwind" of 
the AGB star, was not ejected at very high velocity 
(panel b of Fig. 3). 

4.2. Doppler 6-parameter and synthetic spec- 
tra 

In paper I it was demonstrated that the Red Sys- 
tem (3,0) band is just saturated while the stronger 
(2,0) and (1,0) are saturated. From this it is obvi- 
ous that the Violet System (0-0) band with /v-values 
a factor of ten higher than the Red System must be 
highly saturated. Of course, this circumstance facili- 
tates the detection of the ^"^CN lines. The degree of 
saturation can be judged from Fig. 4 showing curves 
of growth for each N" level. Observed lines from a 
given N" level are combined into a single curve using 
the adopted /-values. The fitted theoretical curve 
is computed from a single line assuming h = 0.51 
km s~^. Clearly the unsaturated lines are all from 
the Red System: 3-0 and 4-0 for the N" = level, 
but extending to the 1-0 band for the N" = 3 level. 
For all relevant A^" levels, the ^^CN/^CN ratio is 
obtainable from a comparison of -"^^CN Violet System 
lines with CN Red System lines of about the same 
strength. Thus, the ^^CN/^'^CN ratio is insensitive 
to the 6-value but dependent on the ratio of the /- 
values for the Violet and Red System. A sensitivity 



to the &-value does exist because, as noted above, the 
precise curve of growth is dependent on the make-up 
of individual CN features. For a detailed analysis we 
must determine the Doppler fe-parametc!r and derive 
column densities from the observed equivalent widths 
by taking into account the details of the blends and 
optical depth effects. 

For each transition (spectral line) there are two un- 
known parameters which affect the observed equiv- 
alent width: the Doppler ^-parameter and the col- 
umn density of the level the transition originates from 
{N{N" j). The ^-parameter must be the same for all 
lines, under the assumption that they are all formed 
from the same gas. The N{N") is the same for those 
transitions for a given level N" of the same isotope. In 
this study we included N" = 0, 1, 2, 3, 4, 5 for CN and 
N" = 0,1,2,3 for "CN, thus ten independent col- 
umn densities. We assumed that the Fi and F2 level 
population is given by their degeneracy level and the 
energy of the average A^" level. We solved for b and 
N{N") in the following manner: 

For a range of 6-parameters (0.10 < 6 < 2.60 km in 
steps of 0.01 km s~^), the column density responsi- 
ble for the observed equivalent width of the blend is 
determined by means of spectrum synthesis. We did 
not fit the observed spectra directly, but instead used 
the observed equivalent width. Although we have ten 
independent measure of the 6-parameter, most of the 
N" levels are only probed by one or two lines and the 
errors are rather large. We have limited the analy- 
sis therefore to the CN TV" = 0,1,2,3 levels which 
are each probed by at least ten blends. For each 
6-parameter, log N{N") and its standard deviation 
where computed. The final 6-parameter is given by 
the combination for which the standard deviation is 
minimal, and by averaging the 6-parameter weighted 
by the number of lines. The CN Doppler-parameter 
was found to be 6 = 0.51 ± 0.04 km s~^. This same 
process was repeated for C2 using the data of Paper I, 
and this resulted in 6 = 0.53 ±0.03 km s~^. Since the 
CN data is of higher quality we adopt 6 = 0.51±0.04 
km s~^ for both CN and C2. Given the 6 value, the 
FWHM of the absorption line profile should be at 
least 0.85 ± 0.05 km s~^. In order to resolve these 
lines a spectral resolution of i? > 350, 000 is needed. 
This confirms that we are not be able to resolve these 
lines profiles in our spectra. The line profiles of the 
synthetic spectra can be perfectly fitted to the ob- 
served ones for this 6-vahie convolved to the spectral 
resolution of the observations. No additional macro- 
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turbulent broadening is needed to explain the obser- 
vations. 

4.3. Column densities and rotational temper- 
atures 

Some additional lines were observed which are 
identified as due to perturbations (Prasad et al. 1992). 
Perturbations occur from the crossing of energy lev- 
els in two different electronic manifolds (Kotlar et al. 
1980). The perturbations in the Violet System (0,0) 
band originate from B^E ?; = levels perturbed by 
the A^n V = 10 level. Perturbations affect the oscil- 
lator strength of lines. Since a study of the pertur- 
bations is beyond the scope of this work, we will not 
include the perturbed lines in our analysis. 

The derived average column density per N" level 
and isotope is listed in Table 4. Given these col- 
umn densities an absolute rotational diagram can 
be constructed (Fig. 5) and the rotational temper- 
ature can be determined. Using a first order fit (con- 
stant Trot), we find Trot(CN) = 12.8 ± 0.6 K, and 
rrot("CN) = 8.0 ±0.6 K. The rotational temperature 
for N" = 0, 1, 2 and 3 of CN is determined very largely 
by the weak CN Red System lines, especially the (3-0) 
and (4-0) lines measured off our McDonald spectra: 
Trot(CN) = 11.1 ± 0.8 K. For N" ^ 4 and 5, the col- 
umn density is set almost exclusively by a CN Violet 
System line which, in the case of N" = 4, is partially 
saturated: the N" = 4 and 5 levels give a temper- 
ature of (rrot(CN) = 10.8 ± 0.8 K) that within the 
errors is identical to the temperature from N" = to 
3. The column density for N" = 4 and 5 seems higher 
by 0.3 dex than given by extrapolation from N" = 
to 3. This offset corresponding to a reduction of the 
Red System /r- values by a factor two (relative to the 
Violet System) is probably excluded by experimental 
and theoretical results on the /-values (see above). A 
single rotational temperature (Table 4) can be fitted 
to the points in Fig. 5. 

Available ^^CN lines, all from the Violet System, 
are partially saturated except for the N" = 3 lines. 
Saturation means that the derived column densities 
are dependent on the b-value and the adopted split- 
ting of the unresolved lines that make up the lines. In 
light of the fact that the 6- values from CN and C2 are 
in excellent agreement, we believe the 6- value is not a 
serious contributor of uncertainty to the determina- 
tion of Trot for ^^CN. Certainly, the column density of 
N" = cannot be lowered (relative to N" = 3) by 0.5 
dex needed to make the ^^CN rotational temperature 



equal to that of CN. 

From a study of millimeter emission lines (e.g. CO 
and CN) of seven planetary nebulae, Bachiller et al. 
(1997b) found a kinetic temperature of the CN line 
forming region of only 25 K, and a rotational temper- 
ature of CN of 5-10 K. Quite unexpectedly we found 
different rotational temperatures for CN and -"^^CN of 
Trot = 11.5 ± 0.6 and 8.0 ± 0.6 K respectively 

In light of the continuing discrepancies between ex- 
perimental and theoretical estimates for the CN Red 
System's /R-values and the radiative lifetimes of the 
A^n state, it is of interest to consider if our spec- 
tra can give useful information on the ratio (/r//v) 
of the Red to Violet System /-values. Of course, the 
fundamental problem is that the Violet strongest lines 
of a given N" are much stronger than the strongest 
Red System lines of the same N" (-"^"^CN lines are of 
comparable strength to the Red System CN lines but 
the ^^CN/^^CN ratio is not independently known). A 
simple check of the /r//v is possible by insisting that 
each band yields the same column density. We have 
conducted this exercise and found that relative to the 
Violet System (0,0) band, the Red System bands give 
a larger total column density by a factor 2.3 (1-0), 2.4 
(2-0), 2.1 (3-0), and 2.0 (4-0). To get consistent col- 
umn densities the /r values have to be decreased by 
this factor. This is clearly unacceptable and we will 
adopt our initial estimated of the /-values to proceed 
with our analysis. It seems that the correction factor 
is wavelength dependent: it increases with increasing 
wavelength. This could suggest that there is scat- 
tered light (contimmm, or CN emission line radiation) 
which fills in the absorption profile. Since our ob- 
servations were collected over several years (the Red 
System was mainly observed in 1992 and the Violet 
System only in 1997), there is the possibility that the 
strength of the lines have decreased with time. This 
idea can easily be tested by simultaneously observing 
both systems. 

4.4. CN/^CN/^CN/Ci^N isotope ratios 

Based on our findings of a different rotational 
temperature for CN and ^^CN, it is clear that the 
CN/^^CN ratio will dependent on the N" levels in- 
cluded in the determination of the column densities. 
The CN/i^CN ratio for each N" level is given in Ta- 
ble 4 with an estimate of its uncertainty. In order 
to determine the isotope ratio of the line forming 
region we add the contributions from all N" levels: 
CN/13CN = E^"=o N{N"f^/ E?,.=o N{N"f^^ = 
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38. Oil. 5. If the measured Tiot is applicable to N" > 5 
levels, these bands contribute insignificantly to the 
summation. 

Upper limits to the equivalent widths of 1.5 mA for 
both the ^^CN and C^^N lines correspond to the fol- 
lowing lower limits: CN/C^^1^{N" = 1) > 2000, and 
CN/"CN(Ar" = 1) > 2000. Assuming that the rota- 
tional temperature of ^^CN and C^^N is equal to the 
one of CN, we find lower limits of CN/C^^N > 2000, 
and CN/i^CN > 2000. 

5. Discussion 

Doppler parameter of 6 = 0.51 ± 0.04 km s~^ sets 
and upper limit on the kinetic temperature of Tkin = 
400 ± 50 K. This is rather high for gas which is about 
10^^ cm from the central star (Meixner et al. 1997). 
Based on radio observations of CN towards post-AGB 
stars, Bachiller et al. (1997b) argue that they find 
Tkin — 10 to 25 K. The fo-parameter must, therefore, 
have an additional contribution from turbulence. The 
CN rotational temperature cannot be used to esti- 
mated the kinetic temperature since the molecule is 
sub-thermal. C2 behaves supra-thermal, but the rota- 
tional temperature derive from the J" = and J" = 1 
level is a good indicator for the gas kinetic temper- 
ature. This would suggest Tkin — 189 K (Bakker & 
Lambert 1997). This value is also much higher than 
would be expected, and it is clear that we need an in- 
dependent accurate determination of the kinetic tem- 
perature. The turbulence required is in the range of 
^^micro = 0-36 to 0.50 km s~^ and depends on the 
adopted Tkin- The source of turbulence is not iden- 
tified, but might be the result of small scale turbu- 
lence, non-Maxwellian velocity distribution, velocity 
stratification, the presence of multiple unresolved ab- 
sorption components due to separate clumps. 

An isotopic ratio as determined from molecules 
('XY and •'XY) does not necessarily refiect the intrin- 
sic isotopic ratio {^X/^ X) of the gas. For example, in- 
terstellar (and probably circumstellar) CO is subject 
to two competing processes which aff'ect this ratio: 
photodissociation by line radiation and ion-molecule 
charge exchange (Watson et al. 1976). The net re- 
sult depends on the local conditions (density, temper- 
ature, and radiation field). 

Since photodissociation for the CO molecule is 
driven by line radiation, the more abundant CO, the 
more the molecule can shield itself against the dis- 
sociating by photons. The CN molecule is photodis- 



sociated by continuum radiation, not by lines. The 
ion-molecule reaction on the other hand will operate 
for CN. The gas-phase ion-molecule reaction for CN 
is given by: 
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C+ + ^^CN ^ 12^+ -I- "CN + AE 



(1) 



where AE is the zero-point energy difference derived 
from the molecular parameters (Prasad et al. 1992) 
and the relation of isotopic dependence of the Dun- 
ham coefficients. We find 31, 58, and 23 K for ^^CN, 
^^CN, and C-^^N respectively. The reaction to the 
right will be referred to as the forward reaction (kf) 
and to the left as the reverse reaction (fc^). C+ can 
be produced by cosmic ray ionization of free carbon 
atoms and molecules containing carbon (e.g. CO), 
and by the reactions of He+ with these molecules 
(where CO is the dominant one). 

Our picture of the chemistry of the CN molecule in 
the circumstellar shell is a simple one. CN molecules 
which are produced by photodissociation of the par- 
ent molecule HCN, are subsequently destroyed by 
photodissociation. Between their formation and de- 
struction, CN molecules participate in the isotopic 
exchange reaction. If the time scale for the isotopic 
exchange reaction is shorter than the lifetime of CN, 
then the exchange reaction will reach equilibrium: 



a = 



with 
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with the square bracket denoting the local density of 
that species in cm~^, R^^ and Rq the equilibrium CN 
and true C isotope ratio. A lower limit to the kinetic 
temperature is given by the CN rotational tempera- 
ture of 12 K. This gives an upper limit on a of 13. For 
higher temperature, a decreases to unity. This clearly 
indicates that if the ion-molecule reaction is in equi- 
librium, we may underestimate the isotope ratio by 
as much as a factor 13! For a more realistic Tki„ ~ 25 
K, and an i?eq = 23, we find a true Rq = C/"C = 80. 

To assess the behavior of the ratio ^^CN/^CN 
in an expanding shell, we use the theoretical work 



11 




Fig. 4. — Curve of growth for each N" level. The displacement of CN (solid line) and ^^CN (dashed line) along 
the horizontal axis gives the isotope ratio derived from that N" level. Squares for the Violet System (0,0), and 
triangles for the Red System bands (triangle pointing up (1,0), right (2,0), down (3,0), and left (4,0)). The over 
plotted theoretical model is for h = 0.51 km s~^. 
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Fig. 5. — Rotational diagram for CN and ^^CN with the column densities of Table 4 (6 = 0.51 km s ^). Note that 
the slope of the two curves are different. This translates in CN having a higher rotational temperature than ^^CN. 
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Table 4: Column densities and the isotope ratios derived from spectrum synthesis for b = 0.51 km s ^, and the 
/-values given in Sect. 3.1 & 3.2. 



AT// 

N 




1 AT/ AT//\ r —2' 

logJy[N ) [cm 








CN/ CJN CN/ °CJN 


CN 











14.80 ±0.17(10) 


13.22 ±0.10(1) 






38.3 ±3.0 




1 


15.05 ±0.17(12) 


13.59 ±0.10(2) 


< 11.72 


< 11.72 


28.8 ±2.1 


> 2000 > 2000 


2 


14.85 ±0.16(17) 


13.16 ±0.10(2) 






48.2 ±2.8 




3 


14.39 ±0.13(11) 


12.44 ±0.10(2) 






89.2 ±3.9 




4 


13.78 ±0.10(01) 


< 11.79 










5 


12.70 ±0.10(01) 












6 


< 11.79 






















pr 


eferred ratios 


Total: 


15.44 ±0.06(52) 


13.86 ±0.09(7) 


< 11.72 


< 11.72 


38.0 ± 1.5 


> 2000 > 2000 


Trot [K]: 


11.5 ±0.6 


8.0 ±0.6 











- numbers in brackets are the number of lines used to determine the average column density. 

- errors are the standard deviations on the average, if less than 4 lines, 0.10 was adopted. 
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Fig. 3. — A comparison of the derived column density 
of the N" = 2 level from the CN Violet System (0,0) 
R(2) line at 3873.37 A. The solid line is obtained from 
the spectrum synthesis method, while the dashed line 
uses the effective oscillator strength and a curve of 
growth for a single, unblended line. The separation 
between the two principal lines of the blend is 0.54 
km s~^. 



of Cherchneff et al. (1993) on the chemistry of 
the circumstellar shell surrounding the carbon star 
IRC +10216. The sheU radius where CN is more 
abundant than half its peak abundance are between 
0.7 X 10^"^ and 1.3 x 10^'^ cm. Together with the ex- 
pansion velocity this gives a lifetime for CN of about 
2 X 10'^ years. Basically this is the lifetime of a single 
CN molecule after photodissociation of HCN, before 
the photodestruction of CN, not taking into account 
any other reaction that could produce or destroy CN. 

The time scale for the ion-molecule reaction is more 
difficult to estimate. Adams et al. (1985) computed 
reaction rates (fc/) for polar molecules and found that 
they are as high as 10^*" cm'^ s~^ for low tempera- 
ture gas (Tkin < 50 K). The radius at which the CN 
abundance peaks (r ~ 1.0 x 10^^ cm^'^) has a C+ 
abundance of 3 x 10~* cm""^. The time scale for the 
isotopic exchange reaction as derived from the for- 
ward rate coefficient and this C^ particle density is 
of the order of 1 x 10^ years. The isotopic exchange 
reaction time scale is comparable to the CN life time. 
This suggest that we have to look into this process in 
some greater detail. 

The effect of the isotopic exchange reaction should 
ideally be studied by computing self consistent chem- 
ical models for the abundances and radiation field ap- 
propriate for HD 56126. Such a study is beyond the 
scope of this work, instead we will make some esti- 
mates. A model of the circumstellar shell of the car- 
bon star IRC +10216 has been computed by Cherch- 
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Fig. 6. — (a) The density of CN and ^^CN (solid lines) and C+ and ^^0+ (dashed lines) as function of distance to 
the star (adopted from Chcrchncff ct al. 1993). (b) The contribution of the column density to the total column 
density, normalized to a peak value of 1. (c) The local isotope ratio for CN and C"^ as function of radius for reaction 
rates of logA;/=-ll to -3 in steps of 1. kf = 1.0 x 10~^ is enhanced, (d) The final observed isotope ratio, with a 
dot at our preferred rate coefficient. 
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nefF et al. (1993). This model uses a large number of 
chemical reactions, the interstellar radiation field, cir- 
cumstellar reddening, and computes the abundances 
of the most important molecular species. The two ma- 
jor differences between IRC +10216 and HD 56126 are 
(i) the inner shell radius of HD 56126 will be much 
larger than for IRC +10216, and (ii) the radiation 
field of HD 56126 (7000 K, Klochkova 1995) is signifi- 
cantly hotter than that of IRC +10216 (3500 K). The 
larger inner radius docs not affect the chemical struc- 
ture of the remaining shell. A 7000 K stellar radiation 
field has fewer UV photons than the interstellar UV 
field (see Paper II) such that the interstellar radiation 
field will dominate the photodissociation processes. 
Therefore the chemical structure of IRC +10216 is 
applicable to HD 56126. From the model, we infer 
that CN is produced mainly by photodissociation of 
HCN and destroyed by photodissociation into C and 
N. From Chcrchneff ct al. we take the CN and C"*" 
density (panel a of Fig. 6). Assuming that the iso- 
topic exchange reaction for CN has reached equilib- 
rium (Tkin = 25 K), we assume a true isotope ra- 
tio of Ro = 67. From this we compute where the 
contribution to the observed column density comes 
from (panel b of Fig. 6): the line is formed where 
the CN density is highest. For a range of reaction 
rates we have computed the local CN and ^^CN den- 
sities, while the C"*" density in unaffected by the iso- 
topic exchange reaction, (panel c of Fig. 6). C+ is 
mainly formed from CO and the processes involved 
in the production of C+ arc not coupled to CN. The 
assumed reaction rate is kf = 1.0 x 10~^. We see 
that the transition from the initial CN isotope ratio 
to the equilibrium value occurs where the CN density 
is highest and that this transition is rather fast. From 
these densities, we can computed the predicted ob- 
served CN isotope ratio (panel d of Fig. 6). To obtain 
an observed CN isotope ratio of 38 at Tki„ = 25 K, we 
have to start with an intrinsic isotope ratio of i?o = 67 
and i?oq = 19. Clearly, the isotopic exchange reaction 
is important. Based on this analysis we suggest that 
the intrinsic carbon isotope ratio is ^^C/^^C=67 for 
Tkin = 25 K. Adams et al. (1985) show that the rate 
coefficient is dependent on the rotational quantum 
number. The higher the rotational quantum num- 
ber, the lower the rate coefficient. If this dependence 
is rather strong, then the isotopic exchange reaction 
would be able to explain the difference in rotational 
temperature between CN and ^^CN: for N" = more 
^^CN is formed than for N" = 5 leading to a lower 



CN isotope ratio as derived from N" = than from 
TV" = 4. 

Isotopic exchange reactions arc very efficient at low 
temperatures and for high dipole moment molecules. 
From this we predict that CO (with a lower dipole 
moment) will be less affected, while C2 with no dipole 
moment is unaffected. C2 will therefore give the true 
isotope ratio and combined with CN this may give the 
kinetic temperature of the gas. The isotopic exchange 
reaction also enhances the abundance of ^^CN and 
C^^N. 

In order to derive the mass-loss rate from the CN 
absorption lines we adopt the CN abundance (in num- 
ber of particles) of Xqn = 1.0 x 10"^ for AFGL 2688 
(Bachiller et al. 1997a). AFGL 2688 (the Egg Neb- 
ula) is so far the only C2CN stars for which ra- 
dio emission of CN has been observed. The inner 
and outer radius of the ejecta of HD 56126 have 
been determined by Meixner et al. (1997) by mod- 
eling the spatial distribution of the infrared source: 
ri„„e, = 45 ± 2 X 10^5 cm (1.0+0.15"), and router = 
155 + 15 X 10^^ cm (3.6 x 3.3" at 11.8 fj,m). Assuming 
that CN coincides with the infrared source, we find 
M = 4.5 + 0.8 X 10"^ Mq yr-i. 

The prototype of the AGB stars, IRC +10216/CW 

Leo, is a massive, highly-evolved carbon star with 
3Mq < MzAMS < 5M0 (Guelin et al. 1995). 
For IRC +10216 the isotope ratios are well deter- 
mined and give an estimate of the ratios which one 
might expect to detect for a carbon-rich post- AGB 
star like HD 56126: G/^^C=Utl, ^'^C/^'^C> 62600, 
N/i-'^N>5300, O/"O=840t?-^^, and 0/1^0=1260^^;^^ 
(see Forestini & Charbonnel for an overview). Our 
estimate of the ^^C/^^C ratio for HD 56126's shell is 
consistent with that of IRC +10216. Results for cir- 
cumstellar shells around four other carbon stars were 
provided by Kahane et al. (1992): C/^^ < qq^ > 39^ 
32tJ." and 3lt^. These too suggest that HD 56126 is 
not exceptional. 

A comparison may also be made with isotopic ra- 
tios pertaining to the photosphere of cool carbon 
(AGB) stars. It must be borne in mind that the 
photospheric ratios may evolve before the stars shed 
the material that becomes the circumstellar shell of 
a post- AGB star. Luminous AGB stars of interme- 
diate mass experience H-burning at the base of their 
convcctive envelope. This decreases the C/^^C ratio 
but rather quickly also converts the star from C-rich 
to O-rich. It seems improbable that C-rich post- AGB 
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stars will have evolved from this intermediate stage. 
At lower luminosity and for lower mass stars, ^^C 
from the He-burning shell is added following a ther- 
mal pulse. Repetition of this addition will increase 
the photosphere's C/^^G ratio unless H-burning con- 
verts ^^C to ^^C. On this scenario, the C/^^C ratio 
of a post-AGB star might be higher than that of a 
typical AGB progenitor. 

Through analysis of infrared spectra of N-type car- 
bon stars, Lambert et al. (1986) showed that most 
stars had a C/^^C ratio in the range 30 to 70: a value 
of 67 is not unusual. Ohnaka & Tsuji (1996), how- 
ever, from analysis of spectra near 8000 Areported, 
systematically lower ratios but a ratio of 67 would 
not be unusual, although on the high end of the dis- 
tribution function. This possible shift between the 
C/^^C ratio of stars and HD 56126 may be due to the 
fact that the stellar ratios continue to evolve before 
the post-AGB phase is reached. On the other hand, 
quantitative analysis of SiC circumstellar grains ex- 
tracted from meteorites shows a C/^^C distribution 
function that mimics well Lambert et al.'s stellar dis- 
tribution, which suggests, perhaps, that evolutionary 
changes are minimal. 

A priority goal has to be the determination of more 
accurate C/^'^C ratios for HD 56126 and an adequate 
sample of post-AGB stars. In the case of HD 56126 
the ratio is presently uncertain for several reasons. A 
minor contributor to uncertainty is the combined use 
of CN Violet and Red System lines and the presence 
of uncertainty over the /-values of the Red System. 
This uncertainty could be avoided by re-observation 
of the strongest CN Red System bands to higher S/N 
ratio such that their ^^CN lines are detectable. A 
key problem will remain in that the CN molecules 
are fractionated in the shell and this has to be mod- 
eled in order to extract the C/^^C ratio that is of 
primary interest. This extraction requires a better 
understanding of physical conditions in the shell. We 
suggest that this understanding may be obtained by 
analyzing spectra of other C-containing molecules. In 
particular, C2 is accessible by the Phillips system, and 
CO is probably detectable in the infrared. Lacking 
an electric dipole moment, the isotopic ex(;hange re- 
action involving C2 will be very slow, and, as we have 
already shown observationally (Bakker & Lambert 
1997), the C2 and CN rotational ladders are influ- 
enced differently by the local density and the ambient 
radiation field. The CO molecule has a smaller dipole 
moment than CN, and its photodissociation is line- 



dominated. Measurements of the ratios C2/^^C^"'C, 
and CO/"CO together with the CN/^^CN ratio will 
shed light on the physical conditions in the shell, and 
on the true ^^C/^'^C ratio. If such data can be as- 
sembled for a sample of post-AGB stars, it will be 
possible to trace a link with their AGB antecedents. 
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A. Honl-London factors and effective oscillator strength 

Honl-London factors were computed from tlic equations given in Schadcc (1964) and Kovacs (1969), and nor- 
malized according to the sum rule ^ j, Sj'j" = (2 — (5o.a'+a") x (25"' + 1) x (2J" + 1) as stated by (Larsson 1983). 
For both the Red and Violet System we take S" — 0.5 since both originate from the same ground level (^S). For 
the Red System (H - S, A' = 1 and A" = 0) we find So,i = and for the Violet System (S - E, A' = and A" = 0) 
•^0,0 = 0. For the Red System we summed over twelve transitions following the selection rule (A J — —1, 0, 1 and 
AN = -2, -1, 0, 1, 2) (six for N" ^ J" + 0.5 and six for N" = J" - 0.5 for a given J") (Pi, Qi, Ri, QPai, ^Q2i, 
PQi2, *^Ri2, P2, Q2, R2, '^R2i and °Pi2) and for the Violet System over six transitions following the selection rule 
(AJ = -1, 0, 1 and AA^ ==-1,1) (three for A^" = J" + 0.5 and three for N" = J" - 0.5 for a given J") (Pi, Ri, 
'^P2i, *^Ri2, P2 and R2). For both the Red and Violet System this results in the relation: J2j' fj'J" = 2 x fv'v"- 
The oscillator strength of individual lines were checked against those of the SCAN tape (for the Red system) and 
those listed in literature (for the Violet system). 

Given a Doppler broadening of b = 0.51 km s~^, transitions which are separated from each other by less than 
0.012 A (Violet System), and 0.022 A (Red System), overlap and the two transitions should be treated as one 
line with an effective oscillator strength. Lines which are separated by more than this interval should be treated 
separately, although they might not be resolved in the spectrum. In physical terminology: if lines are separated by 
less than the Doppler 6-parameters, their optical depth at any frequency point should be added before computing 
the line profile (/ = e"^'^^^'^^^). If lines are separated by more than the Doppler 6-parametcr, the line profiles at 
any frequency point can be added (/ = 0.5 x e~'^^ + 0.5 x e"'^^), although adding the optical depth is also fine. For 
each blend we can determine an effective oscillator strength in the following way: 

The Fl (J" = N" + 0.5) and F2 (J" = N" - 0.5) levels have different degeneracy: 



ge{Fl) ^2x J" + l = 2x N" + 2 (Al) 
ge{F2) = 2x J" + l = 2xN" (A2) 

Each N" levels has therefore a degeneracy of 

geiFl + F2) ^2 X {2x N" + 1) (A3) 

This is even valid for A'^" = although that level has no F2 levels. The reason for this is that each J" level 
is split into 2 x J" + 1 so called m-levels. Those levels are only separated from each other under the influence 
of a strong external magnetic field (Zeeman splitting). In the case of circumstellar CN there is no observational 
evidence to assume that these levels are separated which indicates that there is no measurable magnetic field in the 
line forming region. To first order the m levels for Fl and F2 for a given N" have the same energy and therefore 
the same population. This allows us to write: 



A^(A^", Fl) = ge{Fl)/ge{Fl + F2) x N{N", Fl + F2) (A4) 

N{N", F2) = g^{F2)/g,{Fl + F2) x N{N", Fl + F2) (A5) 

N{N", Fl + F2) = N{N", Fl) + N{N", F2) (A6) 

and 

UN{N") = Y,m{N") (A7) 

i 

h^N{N") = MFl) X (AT" + 1) + E /^(^2) X AT" j x (A8) 
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As examples we take two blends. AN = N' - N" = at 6194.5 A (Red System): 

N(N" = 1) 

fesN{N" = 1) = (0.091 X 2 + 0.363 x 1) x —5^— ^ = 0.182iV(iV" = 1) 

o 

and AN = N' - N" = -1 at 3876.3 A (Violet System): 

N(N" = 2) 

fesN{N" = 2) = (0.0132 X 3 + 0.0022 x 2 + 0.0110 x 2) x ^ '- = 0.0n2N{N" 

o 

B. Extra line lists 



Table 5: Index of refraction for air. Table will not appear in journal. 

Band nair 

Violet System (0,0) 1.000283 

Red System (1,0) 1.000274 

Red System (2,0) 1.000275 

Red System (3,0) 1.000276 

Red System (4,0) 1.000277 
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Table 6: CN Red System (1,0) with /(i,o) = 1-90 x 10"^. Derived column densities are given in Table 4. Table will 
not appear in journal. 



B(N") 


At-ogI- \J\.\ 

' 'I c&l L J 


ffyj' T f\F" J" X 10^ 


feff X WN(N"'\ fcm-21 


Wx ImA] 


Remark 


«R2i(3) 


9129.142c 


1.75 


1.00N(3) 


no 






9133.730c 


1.95 


1.17N(2) 


no 






9136.579c 


2.19 


1.46N(1) 


no 




"R91 CO) 


9139.686c 


2.62 


2.62N(0) 


33.7 




R2(3) 


9141.169c 


3.55 








RQ2i(3) 


9141.189c 


3.77 


3.67N(3) 


35.3 






9141 870c 


3.55 








^091 (2) 


9141.886c 


4.11 


3 89Nf2) 


57.8 




^ \ / 


9142 828c 


3.76 










9142.838c 


4.62 


4.33N(1) 


57.8 


polluted 




9144 042c 


6.34 


6 34N('0'l 


57.8 




09(1) 


9147.208c 


6.34 


3.18N(1) 


61.5 




'^P9l('l) 


9147 217c 


1.60 








09(2) 


9149.167c 


5.53 


3.29N(2) 


47.0 




21 v^/ 


9149 182c 


1.79 








09(3) 


9151.381c 


5.67 


3.45N(3) 


31.1 




QP2i(3) 


9151.401c 


1.79 








P2(2) 


9153.532c 


1.59 


0.636N(2) 


27.4 




P2(3) 


9158.658c 


2.02 


0.866N(3) 






Ri(3) 


9177.032c 


4.86 


2.78N(3) 


31.5 




Ri(2) 


9179.911c 


5.32 


3.19N(2) 


47.6 




Ri(l) 


9183.209c 


6.38 


4.25N(1) 


101.8 


polluted 


Ri(0) 


9186.932c 


10.10 


lO.lON(O) 


59.9 




^Rl2(3) 


9189.465c 


3.72 








Qi(3) 


9189.486c 


5.61 




97.8 


polluted 


^Rl2(2) 


9189.581c 


5.03 








Q2(2) 


9189.597c 


5.14 








^Rl2(l) 


9190.120c 


8.96 


5.83N(1) 


62.9 




Qi(i) 


9190.129c 


4.26 








^Ql2(2) 


9196.511c 


3.38 


1,802N(2) 


42.2 




Pi(2) 


919().525c 


O.To 








^Ql2(3) 


9199.171c 


3.59 


2.26N(3) 


22.2 




Pi (3) 


9199.192c 


1.27 










<J20G.114c 


0..j2 


0.223X(3) 


no 





no: not observed thus not available 
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Table 7: CN Red System (2,0) with /(2,o) = 9.03 x 10"^. Derived column densities are given in Table 4. Table will 
not appear in journal. 



B{N") 


r V . 

Arest [AJ 


fN'J',N"J" X 10'' 


/eff X 10^N{N") [cm-^] 


Wx [mA] 


Remark 


"R2l(0) 


7871.654c 


1.24 


1.24N(0) 


29.4 


polluted 


R2(3) 


7872.889c 


1.68 








^Q2l(3) 


7872.905c 


1.79 


1.74N(3) 


18.5 


polluted 


R2(2) 


7873.332c 


1.68 








^Q2l(2) 


7873.343c 


1.96 


1.85N(2) 


38.8 


polluted 


R2(l) 


7873.985c 


1.79 








^Q2l(l) 


7873.992c 


2.19 


2.06N(1) 


45.8 


polluted 


^Q2l(0) 


7874.852c 


3.01 


3.01N(0) 


39.1 




Q2(l) 


7877.198c 


3.01 


1.51N(1) 


98.2 


polluted 


'^P2l(l) 


7877.205c 


0.76 








Q2(2) 


7878.686c 


2.62 


1.56N(2) 


35.6 




^P2l(2) 


7878.697c 


0.85 








/~\ /'o\ 

Q2(3) 


7880.384c 


2.69 


1.43N(3) 


19.2 




'3P2i(3) 


7880.400c 


0.85 








P2(2) 


7881.889c 


0.76 


0.30N(2) 


nd 




P2(3) 


7885.725c 


0.96 


0.41N(3) 


nd 




Ri(3) 


7899.481c 


2.31 


1.35N(3) 


42.0 


polluted 


Ri(2) 


7901.520c 


2.52 


1.51N(2) 


31.6 




Ri(l) 


7903.892c 


3.03 


2.02N(1) 


41.3 




Ri(0) 


7906.598c 


4.79 


4.79N(0) 


46.3 




^Rl2(3) 


7908.597c 


1.76 








^Rl2(2) 


7908.611c 


2.39 




50.1 


polluted 


Qi(3) 


7908.613c 


2.66 








Qi(2) 


7908.622c 


2.43 








^Rl2(l) 


7908.959c 


4.25 


2.76N(1) 


53.3 


polluted 


Qi(i) 


7908.966c 


2.02 








^Ql2(2) 


7913.692c 


1.60 


0.86N(2) 


19.1 




Pi (2) 


7913.704c 


0.36 








^Qi2(3) 


7915.714c 


1.70 


1.08N(3) 


18.5 




Pi (3) 


7915.729c 


0.61 










7920.8()4c 


0.2.J 


0.11X(3j 


11(1 





no: not observed thus not available 
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Table 8: CN Red System (3,0) with /(3,o) = 3.34 x 10"^. Derived column densities are given in Table 4. Table will 
not appear in journal. 



Tit AT"\ 


Arost [AJ 


JN'J',N"J" X lU 


Jeff- X lU iv (iv ) [cm J 


Wx [mAJ 


Remark 


'^■i?2i(4) 


6916.8681 


0.282 


0.157N(4) 


no 








U.310 


A 1 '7'71VT ^ 0\ 

L).i771N(^5j 


no 








r\ AO 
0.343 


(J.zUblN [2) 


no 




it2l(l) 


6922.5681 


0.387 


0.258N(lj 


no 






0924. OZ7i 


U.4d1 


L).4bilN ({}) 


1 A T _L 1 A 
10.7 ± 1.0 




i?2(4) 


6925.811c 


0.627 








<y2i(4) 


t^noc 001 1 

byzo.ozll 


0.612 


A 1 A AT ^ A\ 

U.biyiN (4j 


1 n 1 A rr 
1.9 ± 0.5 


tentative 


^?2(3) 


6925.893c 


0.619 












U.DDO 


U.044iN yo ) 


0.0 ± U.5 




n2[2) 


DyzD.iy4c 


0.622 








Q2l(2) 


6926.1781 


0.725 


0.684N(2) 


16.3 ± 0.5 




il2(lj 


byZD-DODC 


U.boo 








<y2i(l) 


byZD.DOUl 


U.oii 


0./60iN(lj 


00 A _l_ A C 

z3.0 ± 0.5 




y2i(Uj 


Dyz7.c)U41 


1 1 1 A 

i.iiU 


111 AM/ rt\ 

I.IIOIN (Oj 


AO C\ _L C A 
4<3.9 ± 5.0 


polluted 


/~i / 1 ^ 

^2(1) 


byzy.ubzi 


1 1 1 A 


0.555iN(lj 


100 _L A C 

IS. 8 ± 0.5 




^-r2l(l) 


69/9. 0o7c 


A A'VO 

0.278 








^2(2) 


693U.Z / 21 


0.969 


0.5761X1 (zj 


1 /I A _L A C 

14. U ± 0.5 




^-r2l(/) 


693U.oUzc 


A 1 /I 

0.314 










6931. OMl 


A AAO 

0.992 


A /? A C AT / \ 

0.605iN (3j 


a i\ 1 1 A 

6.U ± 1.0 


line too weak 


^^^21 (3) 


G931.663C 


0.315 








-^2(2) 


6932.7151 


0.278 


0.111N(2) 


6.0 ±0.5 




<32(4) 


6933.1791 


1.030 


0.627N(4) 


< 1.5 




^^21 (4) 


6933.206c 


A A c: 

0.305 








-^2(3) 


6935.7201 


0.353 


0.151N(3) 


< 1.5 




i^(4j 


6938.9101 


OAA 

0.399 


A 1 ^^AT/ 1 \ 

0.1 / /i\(4j 


^ 1 c 
< 1.5 




i?l(4) 


6945.2731 


0.811 


0.451N(4) 


< 1.5 






by4b.4 /yi 


A r A 

0.850 


A A 0£i AT ^ 0\ 


C A _L 1 A 
5.9 ± 1.0 




i?l(2) 


6947.9621 


0.929 


0.557N(2) 


13.0 ± 0.5 


polluted 




6949. m31 


1.115 


r\ 'T AO AT / 1 \ 

U.74cilN(ij 


/I A _l_ 1 A 

z4.U ± l.L) 




i?l(0) 


6951.8161 


1.766 


1.766JN(0j 


24.0 ± 5.0 




^^12(2) 


6953. 4i5c 


A 001 

0.881 










6953.4611 


0.898 




/I 1 1 AC 

24.1 ± 0.5 






Dyoo.4oyc 


U.DOO 








Vu3j 


cnco ACM 

6953.4011 


A noA 
0.9(50 










6953 643c 


1.568 


1.021N(1) 


25 7± 5 




Qi(i) 


6953.6431 


0.747 








'^i?i2(4) 


6953.776c 


0.532 








Qi(4) 


6953.7871 


1.040 


0.814N(4) 


2.3 ±0.5 


tentative 


'^<3l2(2) 


6957.302c 


0.593 


0.316N(2) 


8.4 ±0.5 




-Pi (2) 


6957.3011 


0.132 








^Ol2(3) 


6958.905c 


0.631 


0.397N(3) 


4.6 ±0.5 




A (3) 


6958.9071 


0.222 








^Ql2(4) 


6960.782c 


0.612 


0.433N(4) 


< 1.5 




A (4) 


6960.8051 


0.289 








^Pl2(3) 


6962.7951 


0.092 


0.039N(3) 


< 1.5 




^Pl2(4) 


6966.2121 


0.138 


0.061N(4) 


< 1.5 





no: not observed thus not available 



24 



